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ABSTRACT 

We present low resolution UV-blue spectroscopic observations of red giant stars in the globular 
cluster M53 ([Fe/H]= —1.84), obtained to study primordial abundance variations and deep mixing 
via the CN and CH absorption bands. The metallicity of M53 makes it an attractive target: a bimodal 
distribution of 3883 A CN bandstrength is common in moderate- and high-metallicity globular clusters 
([Fe/H] > -1.6) but unusual in those of lower metallicity ([Fe/H] < -2.0). We find that M53 is an 
intermediate case, and has a broad but not strongly bimodal distribution of CN bandstrength, with 
CN and CH bandstrengths anticorrelated in the less-evolved stars. Like many other globular clusters, 
M53 also exhibits a general decline in CH bandstrength and [C /Fe] abundance with rising luminosity 
on the red giant branch. 

Subject headings: star clusters and associations 



1. INTRODUCTION 

Star-to-star variations of the strength of the 3883 A CN 
band are nearly universal in moderate-metallicity ([Fe/H] 
> —1.6) globular clusters of the Milky Way, but unusual 
in low-metallicity globular clusters ([Fe/H] < —2.0) (see 
Gratton, Sneden & Carretta 2004 for a thorough review). 
Typically a CN bandstrength is measured as a magni- 
tude difference between the flux within the 3883 A CN 
band and the flux in a nearby comparison band chosen 
to be free of CN absorption. The index 5(3839) intro- 
duced by Norris et al. (1981) is one such measure of CN 
band strength. Not only does 5(3839) vary from star to 
star within individual globular clusters, its distribution 
is often bimodal in moderate-metallicity ([Fe/H] > —1.6) 
clusters. It is often anticorrelated with the strength of 
CH absorption in the G band at 4300 A among giants 
of similar absolute magnitude. This anticorrelation in- 
dicates that CN-strong stars are depleted in carbon. In 
order to have strong CN bands with low [C/Fe] abun- 
dances, the CN-strong stars must also be enhanced in 
nitrogen. This abundance pattern has directed inves- 
tigations into the origin of the bimodality toward stel- 
lar sources, where hydrogen shell burning can produce 
abundance enhancements and depletions such as those 
observed in globular clusters in C and N (and also O, 
Na, Mg, and Al). Meanwhile, few CN-strong giants are 
known in the halo field (see, e.g., Pilachowski et al. 1996; 
Shetrone et al. 1999), indicating that the enrichment 
process responsible for CN bimodality only operated in 
cluster (or proto-cluster) environments. This paper con- 
siders whether the apparent lack of CN bimodality in 
low-metallicity globular clusters is a true lack of abun- 
dance variation or simply an effect of low overall metal 
abundance. 

Previous investigations into the behavior of CN abun- 
dance at low metallicity have been inconclusive as to 
whether large variations in bandstrength can be present. 
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Trefzger et al. (1983) and Shetrone et al. (1999) re- 
port several CN-strong red giants with large nitrogen 
abundances in M15 ([Fe/H]= —2.26), and Carbon et 
al. (1982) found variations of [C/Fe] and [N/Fe] but no 
strong 3883 A CN bands among the red giant stars of 
M92 ([Fe/H]= -2.28). Briley et al. (1994) showed that 
the carbon and nitrogen abundances measured within 
M92 would result in a bimodal distribution of 5(3839) 
if the overall metallicity were scaled up to [Fe/H]= —1.0 
dex. Smith & Norris (1982) and Briley et al. (1993) both 
studied M55 ([Fe/H]= -1.81, Harris 1996), and found 
that it is quite homogeneous in CN bandstrength. That 
lack of CN-strong stars could be the result of low over- 
all metallicity inhibiting the formation of CN molecules, 
low [C/Fe] among the brighter red giants, or a sign that 
whatever enrichment process caused CN variations in 
moderate-metallicity globular clusters did not operate at 
lower metallicities. A study of additional globular clus- 
ters with metallicities similar to M55 could help decide 
between these alternatives. We chose M53 as the subject 
of this study because it is similar in metallicity to M55 
(based on our isochrone fits discussed below), and will 
hopefully provide insight into whether low-metallicity 
globular clusters contain CN-enhanced stars. While our 
study does not comment directly on the process that 
causes abundance inhomogeneities in globular clusters, 
we note here that a lack of abundance variation in the 
lowest-metallicity clusters would provide a constraint on 
all hypotheses about that process. 

To understand the CN and CH index behavior, we 
compare spectroscopy of the 3883 A CN band for red 
giants in M53 to the Norris et al. (1981) study of NGC 
6752, a moderate-metallicity globular cluster ([Fe/H] = 
— 1.56) that shows clear CN bimodality. We also convert 
CH bandstrengths into [C/Fe] abundances using syn- 
thetic spectra and look at the relations between carbon 
abundance and CN bandstrengths. An intrinsic spread 
is found in CN bandstrength, but it is smaller than the 
range observed in NGC 6752. 

This result is interesting in the light of work by 
D'Antona et al. (2002), who use stellar evolution models 
to draw connections between helium abundance and hor- 
izontal branch morphology. They find that globular clus- 
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Fig. 1. — Color-magnitude diagram for the red giant branch of 
M53. The photometry is from Rey et al. (1998) for stars identified 
as red giants in Cuffey (1965). Filled circles are for stars observed 
with Kast, open circles are stars we did not observe, crosses are for 
stars dismissed in §2 as nonmembers. 

ters with solar helium abundances have short, red hori- 
zontal branches, while extended blue horizontal branches 
result from the smaller core mass of helium-enhanced 
stars. In their model, the helium enrichment comes from 
AGB feedback, implying that globular clusters with red 
horizontal branches, like M53, should exhibit only minor 
CNO-product enrichment. 

2. THE DATA SET 

Our sample is comprised of 24 red giant stars in M53, 
selected from the photometric atlas of Cuffey (1965). In 
order to have a data set as homogeneous and intercom- 
parable as possible, all of the observations were made 
using the Kast double-beam spectograph on the Shane 
3-meter telescope at Lick Observatory. A mirror was 
used in place of a dichroic to direct all the stellar light 
to the blue side of the spectrograph, and the 452 1/mm 
grism and a slit width of 1.5" produced a mean disper- 
sion of 2.54 A /pixel and a mean resolution of 6.3 A. 
Exposure times varied from 1800 to 5400 seconds, giving 
a signal-to-noise ratio per pixel at the 3883 A CN band 
between 20 and 60. 

Photometry for the program stars was taken from the 
CCD study of Rey et al. (1998). Cross-identification 
was done by taking the full online data set 4 , plotting the 
locations of all stars brighter than 16th magnitude, and 
visually comparing that plot with the Digitized Sky Sur- 
vey i?-band image of M53. Figure 1 is a color-magnitude 
diagram for the stars observed in M53, showing the rea- 
sonably uniform coverage of our data set over a range in 
absolute magnitude of —0.5 < My < —2.5. We adopt an 
apparent distance modulus of (m — M)y — 16.31 from 
the online catalog of Harris (1996). All but four of the 
stars observed are consistent photometrically, spectro- 
scopically, and in radial velocity with being first-ascent 
red giant branch stars. M53 I, SI R6 5, S3 R2 1, and S3 
R2 15 (shown as crosses in Figure 1) are excluded from 
the sample because it is clear from their radial velocities 
and spectra that they are not cluster members. 

Table 1 lists identifiers, positions, dates of observation, 
photometry, exposure times, three bandstrength indices, 
and measured carbon abundances for each star in our 
sample. The X and Y offsets refer to the distance, in 
arcseconds, of each star from cluster center, and are given 
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Fig. 2. — Nightly offsets between CH and CN bandstrength 
indices measured from extinction-corrected, unfluxed Kast spectra 
of flux standard stars and those measured from HST spectra (taken 
from the CALSPEC database at STScI, Bohlin 1996) of the same 
flux standards. On 2001/04/15, the standard star observed was 
Kopff 27; on all other nights the standard star was Feige 34. 

in the online data associated with Rey et al. (1998). 
The results from the spectroscopy are listed as 5(3839), 
S(CH), and (777(4300). The index 5(3839) is a measure 
of the 3883 A CN band, and S(CH) and C7J(4300) (the 
later of which was introduced by Harbeck et al. 2003) 
both represent the depth of the 4300 A CH feature. We 
define S{CH) as 



S(CH) = -2.5 log - 



f 4320 T 
J 4280 



J4050 ^ dX + J4330 J ^ dX 

where I\ is the number of counts per pixel at a given 
wavelength. Values of CH(4300) are given to facilitate 
comparisons with previous work. 

The spectroscopic data were reduced using the XIDL 
suite of IDL programs developed by Jason Prochaska and 
Gabriel Prochter. 5 The Kast-specific programs were used 
to identify observation types, create a nightly HeHgCd 
composite arc spectrum and dome flat-field image, do 
flat-field corrections, trace, extract and wavelength- 
calibrate one-dimensional spectra, and combine multiple 
spectra of the same object. The XIDL routines did not 
include an option for correcting spectra for atmospheric 
extinction; since some of our spectra were taken at an 
airmass of 2, a script was written to perform those cor- 
rections prior to the measurement of indices from the 
data. The correction was based on the mean extinction 
curve given on the Lick Observatory web site. 

We transformed the bandstrength indices onto a stan- 
dard system defined by indices derived from HST spectra 
of flux standard stars taken from the STScI CALSPEC 
database described by Bohlin (1996). These "standard" 
index values were compared to indices measured from 
the extinction-corrected (but unfluxed) Kast spectra that 
were obtained on each night of one or more of the CAL- 
SPEC flux standards. Figure 2 shows the differences 
(in the sense Shst — Skast) f° r the indices 5(3839) and 
S(CH) for each night of observation. They are quite con- 
sistent from night to night, both within an observing run 
and between runs separated months apart. These diffcr- 



4 Available at http:/ /vizier. u-strasbg.fr/viz-bin/VizieR?-source= J/AJ/lfcAl#6i£i/blc from http:// www.ucolick.org/ xavier/IDL/ 
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TABLE 1 



Star ID a 


X offset (")" 


Y offset (") B 


Date 


T 7 Vl 

V a 


(B — V) b 


t (s) 


5(3839) 


S(CH) 


Cff(4300) 


[C/Fc] 


K 


-139.7 


-85.24 


041 52001 


13.78 


1.50 


2400 


0.1220 


0.6081 


0.5199 


-0.35 


S4 R4 16 


56.59 


-269.46 


05162001 


13.95 


1.30 


3600 


0260 


5693 


0.4858 


-0.65 


G 


-261.67 


99.70 


041 52001 


14.08 


1.32 


4200 


0.0927 


638^ 


0.5412 


-0.35 


S4 R4 24 


-96.74 


-252.57 


05072002 


14.34 


1.20 


3600 


-0.0403 


6639 


0.6180 


-0.35 


S3 R2 22 


-114.81 


-4.25 


05082002 


14.34 


1.17 


1800 


0.0961 


0.5706 


^580 


-0.80 


S4 R5 1 


316.48 


-281.80 


051 62001 


14.56 


1.08 


3600 


-0 0868 


6693 


0.6719 


-0.40 


M 


-102.07 


-163.22 


041 52001 


14.59 


1.08 


2400 


0.0198 


0.6941 


0.7028 


-0.30 


f 


144.35 


-285.43 


051 62001 


14.62 


1.07 


3600 


-0.0827 


0.7090 


69^5 


-0.20 


z 


23.83 


242.97 


05082002 


14.73 


1.08 


3600 


0.0831 


6680 


0.6421 


-0.40 


N 


6.19 


-187.29 


05072002 


14.74 


1.05 


3600 


-0.0664 


0.7417 


0.7363 


-0.10 


S3 R4 23 


-278.38 


151.07 


041 52001 


14.77 


1.09 


2400 


-0.0313 


0.7194 


0.7081 


-0.20 


p 


104.23 


-197.34 


051 62001 


14.83 


1.03 


3600 


-0.0758 


6699 


0.6379 


-0.45 


S4 R4 25 


-104.09 


-2^2 82 


05072002 


14.85 


1.01 


4500 


-0.0292 


0.6401 


^963 


-0.55 


V 


-109.99 


210.18 


05092002 


14.88 


1.04 


3600 


-0.1208 


0.7124 


0.7076 


-0.25 


S3 R2 3 


-77.70 


-105.43 


05082002 


15.08 


0.97 


3600 


-0.1037 


0.6821 


0.6835 


-0.40 


Y 


-8.79 


273.13 


05082002 


15.08 


0.96 


3600 


0.0454 


0.6996 


0.6952 


-0.35 


S4 R2 15 


19.43 


-150.88 


05072002 


15.21 


0.95 


3600 


-0.0142 


0.7138 


0.7645 


-0.25 


SI R2 3 


125.25 


71.33 


05082002 


15.24 


0.97 


3600 


-0.0867 


0.7252 


0.7106 


-0.25 


S3 R5 4 


-290.14 


-264.51 


05092002 


15.25 


0.98 


3600 


-0.1144 


0.7983 


0.8104 


0.10 


S4 R2 4 


71.17 


-131.34 


05092002 


15.30 


0.94 


3600 


-0.0315 


0.7085 


0.7457 


-0.30 


S4 R4 22 


-71.44 


-276.28 


05092002 


15.31 


0.93 


3600 


-0.0418 


0.7017 


0.7473 


-0.35 


S4 R4 27 


-113.07 


-277.44 


05092002 


15.42 


0.93 


1800 


-0.1319 


0.7170 


0.7409 


-0.25 


S4 R4 26 


-112.30 


-267.91 


05072002 


15.47 


0.93 


5400 


-0.1237 


0.7603 


0.7621 


-0.10 


O 


55.48 


-192.49 


05072002 


15.69 


0.89 


3600 


-0.1223 


0.7725 


0.8202 


0.00 



a Identifier from Cuffey (1965) b From Rey et al. (1998) 

ences were applied as offsets to correct the Kast indices 
from the M53 spectra onto a standard system. 

One-sigma errors in the index measurements were esti- 
mated by measuring 5(3839) and S(CH) from individual 
spectra of stars observed more than once. Typically each 
observation of a star consisted of two or three consecu- 
tive CCD exposures. In cases where three observations 
were obtained of each star, we calculated crg( 383 g) as the 
standard deviation in the mean of the individual index 
values, and os(ch) analogously. In cases where only two 
observations were made per star, 0s(3839) was calculated 



01 

to 



as 



0.89xAS(3839) 



, where AS'(3839) refers to the difference 

between the two available values of the CN index, and 
(7 S(CH) analogously. The resulting mean one-sigma index 
errors are (cr s(3S39) ) = 0.02 and (cr s(CH )} = 0.007. 

Figure 3 shows 5(3839) versus My for our data set 
(top panel) and NGC 6752 (bottom panel). NGC 6752 
has a metallicity of —1.56 (Harris 1996), and its 5(3839) 
distribution is typical of moderate-metallicity globular 
clusters, with two distinct sequences separated by a gap. 
The data for NGC 6752 come from Norris et al. (1981). 
The increase in bandstrength as temperature decreases is 
partly caused by more abundant molecule formation at 
lower temperature, but is made steeper by the changing 
shape of the stellar continuum. The M53 data include 
the error bars calculated in the previous paragraph, and 
it is clear that the range in 5(3839) at a given My is 
significant relative to 173(3839) 1 indicating that the width 
of the CN variation is intrinsic to the stars rather than a 
result of measurement errors. We identify two groups of 
M53 red giants in Figure 3 on the basis of their 5(3839) 
indices: a CN-weak group (plotted as open circles) and 
a CN-enhanced group (filled circles). 

3. THE MODELS 

Conversion from the (V, B — V) plane to the (T c ff, 
log g) plane was accomplished by comparing the M53 
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Fig. 3.— The 3883 A CN bandstrength index 5(3839) versus 
absolute My magnitude for red giant stars in M53 (top panel) 
and NGC 6752 (bottom panel). The upward linear trend exists 
because the shape of the stellar continuum and the abundance of 
CN molecules change with temperature. Relatively CN-enhanced 
stars are shown as filled circles. The straight line in the top panel 
is the baseline against which the quantity 55(3839) is measured. 



fiducial sequence from Rey et al. (1998) to isochrones 
from the Victoria-Regina models (VandenBerg et al. 
2006). The 12 Gyr, [Fe/H]= -1.84, [a/Fe]=+0.3 
isochrone was the best match. MARCS model atmo- 
spheres (Gustafsson et al. 1975) were created for 27 (T e g, 
log <?) pairs between My = — 0.5 and My — —2.6 on the 
red giant branch of that isochrone. The model parame- 
ters were set for consistency with the model spectra used 
in Briley & Cohen (2001): [0/Fe]= +0.2, 12 C/ 13 C = 10, 
and Vturb — 2 km/s. The SSG spectrum synthesis code 
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Fig. 4. — Synthetic spectra in the region of the 4300 A 
CH G band. The spectrum of M53 star G is also shown, ver- 
tically offset from the models, for comparison. Parameters are 
[C/Fe]= (-1.0,-0.5,0.0), [N/Fe]= +0.6, [0/Fe]= +0.2, My = 
— 1.5, 12 C/ 13 C = 10, v tur i, = 2 km/s. Vertical lines show the loca- 
tions of the two sidebands (4050 A to 4100 A and 4330 A to 4350 
A ) and the CH absorption bandpass (4280 A to 4320 A ) of the 
S(CH) index. 
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Fig. 5. — Synthetic spectra in the region of the 3883 A CN band. 
The spectrum of M53 star G is also shown, vertically offset from 
the models, for comparison. Parameters are [N/Fe]= (0.0, 0.5, 1.0), 
[C/Fe]= —0.3, and all others the same as in Figure 4. Vertical lines 
show the locations of the comparison sideband (3883 A to 3916 A 
) and the CN absorption bandpass (3846 A to 3883 A ). 

(Bell, Paltoglou, & Trippico 1994 and references therein) 
was used to create model spectra for each of the 27 (T c ff , 
logg) points with [C/Fe] varying from —1.0 to +0.5 dex 
and [N/Fe] varying from —0.5 to +1.0 in steps of 0.1 
dex. Spectra were smoothed to a resolution of 5.4 A and 
a pixel spacing of 1.8 A to better match the data. 

Figure 4 shows synthetic spectra which illustrate the 
strong dependence of the 4300 A G band on carbon abun- 
dance at My = -1.5 (with a fixed [N/Fe]). The verti- 
cal dotted lines delineate the sidebands while solid lines 
define the science band for the S(CH) index. For com- 
parison, the spectrum of M53 G is also shown, vertically 
offset from the synthetic spectra. We chose the sidebands 
to lie close to the science band, but with as little CH ab- 
sorption as possible. The red sideband was necessarily 
narrow, as there is very little [C/Fe]-independent spec- 
trum redward of 4320 A. The blue sideband is farther 
from the science band than we would prefer, but there 
are complex CN molecular features near 4200 A. Posi- 
tioning of the blue sideband (4050 A to 4100 A ) is a 
compromise between band width, location, and [C/Fe] 
dependence. 

Figure 5 shows the dependence of the 3883 A CN band 
on nitrogen abundance at My = —1.5 (with a fixed 
[C/Fe]), and the wavelength bands for the 5(3839) in- 
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Fig. 6. — Generalized histogram of <5S(3839) for M53 red giant 
stars (solid line). The dotted line is the analogous curve for NGC 
6752 red giants (data from Norris et al. 1981). Each histogram 
is normalized to a maximum of unity. A generalized histogram is 
the sum of unit Gaussian curves centered at each data point, with 
the dispersion a of each Gaussian equal to the error on the points 
(which was taken to be 0.02 for M53 and 0.045 for NGC 6752). 

dex are marked with vertical lines. Only the red edge of 
the comparison band is shown as a dotted line because 
its blue edge coincides with the red limit of the science 
band. For comparison, the spectrum of M53 G is also 
plotted, vertically offset from the three synthetic spec- 
tra. The red sideband for 5(3839) is far easier to place 
than either of the sidebands for S(CH), as there is a 
wide region redward of the 3883 A bandhead that is in- 
sensitive to [N/Fe]. The strong Ca II K line is avoided. 
It is impractical to place a comparison sideband to the 
blue of the 3883 A CN feature, since the intrinsic flux of 
a red giant drops precipitously into the UV. There is also 
a forest of metal absorption lines and hydrogen B aimer 
lines in the UV, and furthermore it is problematic to 
find a wide spectral region to the blue that is reasonably 
independent of [C/Fe] and [N/Fe]. The single sideband 
is partially responsible for the strong temperature trend 
in 5(3839), and can also cause significant instrument-to- 
instrument offsets. For these reasons Norris et al. (1981) 
introduced the quantity 55(3839), defined as the mea- 
sured 5(3839) minus a baseline value that is determined 
from a linear fit to the lower envelope of the CN-weak 
stars in the (5(3839), My) plane. The baseline adopted 
for the M53 red giants is 5 = -0.191 - 0.080M y . 

4. ANALYSIS 

We are interested in comparing the CN and CH in- 
dex behavior in M53 to that in globular clusters with 
more moderate metallicity, and use two separate meth- 
ods for making this comparison. In an index analysis, 
we compare the behavior of the 5(3839) and S(CH) in- 
dices in M53 with the variations observed in NGC 6752, a 
typical example of a CN-bimodal globular cluster. In an 
abundance analysis, we determine [C/Fe] abundances for 
individual M53 red giants by comparing data to isoabun- 
dance loci in the S(CH) versus My plane. 

4.1. Index analysis 

NGC 6752 is an excellent example of bimodal CN 
bandstrength and anticorrelated [C/Fe] and [N/Fe] be- 
havior. Figure 6 shows a generalized histogram of 
(55(3839) for M53 red giants (solid line) with the analo- 
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Fig. 7.— The 4300 A CH bandstrength index S(CH) versus 
My, with CN-enhanced stars (identified in Figure 3) plotted as 
filled circles. S(CH) declines with decreasing My, and at magni- 
tudes of My > —1.4 the CN-enhanced stars are slightly offset as 
a group in S(CH) from CN-weak stars. The error bars on indi- 
vidual data points correspond to the measurement errors cs(CH) 
as described in §2. The large crosses represent the mean values 
and standard deviations for the data in three luminosity bins: 
-0.5 > My > -1.2, -1.2 > My > -1.9, and -1.9 > My > -2.6. 
Isoabundance contours from synthetic spectra are shown as solid 
lines. Each curve has the parameters [N/Fe]= +0.6, [0/Fe]= +0.2, 
-2.6 < My < -0.5, 12 C/ 13 C = 10, and v tU rb = 2 km/s. [C/Fc] 
varies from —1.0 dex (bottom line) to +0.5 dex (top line), with a 
spacing of 0.1 dex. It can be seen from the models that S(CH) is 
sensitive to temperature in stars brighter than My ~ —1.4. Mean 
S(CH) values for the three luminosity bins decline faster than can 
be attributed to temperature effects alone, indicating that carbon 
is progressively depleted as stars ascend the giant branch. 

gous NGC 6752 data from Norris et al. (1981) overplot- 
ted as a dotted line. The solid histogram is a sum of 
Gaussian curves with a a of 0.02, the standard error on 
5(3839) for the Kast M53 spectra, while the Gaussian a 
used by Norris et al. (1981) is 0.045. The distribution of 
CN bandstrength in M53 is clearly not well-described by 
a single symmetric Gaussian curve; indeed, a two-sided 
Kolmogorov-Smirnov test comparing (55(3839) for stars 
identified as CN-weak and CN-enhanced in Figure 3 re- 
turns a probability of 2.67 x 10~ 5 that the two groups 
were drawn from the same underlying distribution. How- 
ever, there is a clear difference between the M53 CN dis- 
tribution and that of NGC 6752: the latter has a large 
separation between two clearly distinct CN groups, and 
also has a more well-populated CN-strong group than 
the former. 

Figure 7 shows S(CH) versus My for our M53 data, 
with stars identified as CN-enhanced (see Figure 3) plot- 
ted as filled circles. There is a general decline in S{CH) 
with My, and there is also an anticorrelation between 
CN bandstrength and S(CH) among the fainter stars 
(My > — 1.4). In the later case, the CN-enhanced gi- 
ants have values of S(CH) at the low end of those found 
among the CN-weak giants, although there are only four 
CN-enhanced giants in our sample with My > —1.4. 
Among brighter giants with My < —1.6 there is no 
such trend. The large crosses represent the mean val- 
ues and standard deviations of the M53 data binned 
into three luminosity groups: —0.5 > My > —1.2, 
-1.2 > My > -1.9, and -1.9 > My > -2.6. We have 
overplotted model isoabundance lines for [C/Fe] vary- 
ing from —1.0 dex (bottom) to +0.5 dex (top). All of 
the models plotted here were computed for the param- 
eters [N/Fe]= +0.6, [0/Fc]= +0.2, 12 C/ 13 C = 10, and 
Vturb — 2 km/s. The comparison between the data and 
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Fig. 8.— Left panel: The 3883 A CN bandstrength index 
55(3839) versus S(CH) for all red giants in the sample, with stars 
identified as CN-enhanced in Figure 3 shown as filled circles. The 
three crosses represent the mean and standard deviation in the 
same luminosity bins introduced in Figure 7. There is evidence 
for some deep mixing: stars in the brightest bin, on average, have 
smaller S(CH) and higher 8S(3839) than stars in the two fainter 
bins. Right panel: <5S(3839) versus S(CH) for red giants with 

My > -1.4. 



the isoabundance curves reveals two main features. First, 
at a fixed [C/Fe], the model curves show a clear tempera- 
ture effect on S(CH) for stars brighter than My = — 1.4 
(T GfT < 4500 K). Second, the decline in S{CH) with ris- 
ing luminosity in the data is steeper than the tempera- 
ture effect alone, and this implies that there is a decline 
in [C/Fe] with decreasing My among the red giants in 
M53. 

Figure 8 shows <55(3839) versus S(CH) for our M53 
data. The left plot shows the full data set, and the 
right panel shows only the stars fainter than My = 
— 1.4. There is no clear anticorrelation in the left panel, 
whereas the right panel shows that among stars with 
My > —1.4 those with enhanced CN bands have rela- 
tively low S{CH) indices, as noted also from Figure 7. 
The lack of a more prominent anticorrelation between 
CN and CH in the full data set may be partly due to 
a reduced sensitivity of the index 5(3839) to nitrogen 
abundance at the low metallicity of M53. This point is 
illustrated by Figure 9, which shows the 5(3839) mea- 
surements for M53 versus My, with model isoabundance 
curves for [C/Fe]= —0.3 dex and [N/Fe] ranging from 
—0.5 to +1.0 dex overplotted. Since (as we show below) 
not all of the stars have [C/Fe]= —0.3 dex, the distribu- 
tion of data points in this plot is not necessarily indica- 
tive of true nitrogen abundance. At the location of our 
faintest star in the (My, 5(3839)) plane, changing [N/Fe] 
by 0.1 dex changes 5(3839) by 0.02. Since a factor of two 
change in nitrogen abundance only produces a few-sigma 
change in 5(3839), its usefulness as an indicator of [N/Fe] 
is limited at the metallicity of M53. 

Since 5(3839) has a diminished sensitivity to [N/Fe] 
at the metallicity of M53, we do not attempt to use it 
to measure [N/Fe] for individual stars. However, we can 
comment on the range of [N/Fe] in our sample, and com- 
pare it to the range reported by Norris et al. (1981) be- 
tween red giants in NGC 6752. Figure 10 shows 5(3839) 
versus My for the Kast M53 data, with two sets of mod- 
els overplotted: the lower line is intended to represent a 
typical CN-weak star with [C/Fe]= -0.3 and [N/Fe]= 0. 
This abundance combination was chosen to produce a 
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Fig. 9.— The 3883 A CN bandstrength index 5(3839) versus My 
for stars in M53. Error bars are the measurement errors o"s(3839) 
as described in §2. Isoabundance curves are shown as solid lines; 
their parameters are [C/Fe]= -0.3, -0.5 < [N/Fe] < 1.0, and all 
others the same as in Figure 7. This plot is meant to illustrate the 
dependence of 5(3839) on nitrogen abundance at low metallicity: 
according to the models, there is only a small change in 5(3839) 
for a given change in [N/Fc] particularly in the region of this figure 
occupied by the CN-weakcst giants. Since the model isoabundance 
contours plotted here have [C/Fe] set to —0.3 dex, and not all of 
the M53 red giants have such a carbon abundance, this plot is 
not meant to be indicative of [N/Fe] abundance for any star in 
particular. 

model locus that follows the lower envelope of the M53 
stars in Figure 10. The upper line represents typical 
CN-strong stars with [C/Fe]= -0.6 and [N/Fe]= 0.9. 
The abundance differences between the CN-weak and 
CN-strong loci were chosen to match those reported in 
NCG 6752 by Norris et al. (1981). They measured 
the abundance differences between a pair of CN-weak 
and CN-strong giants in NGC 6752 with My — —1.3 as 
A[N/Fe]=0.9, A [C/Fe] = -0.3. Similar differences were 
found by Da Costa & Cottrell (1980) for another pair of 
giants in NGC 6752 that are one magnitude fainter. The 
two isoabundance curves in Figure 10 bound most of the 
data points for M53, meaning that while some pairs of 
stars in this cluster may have carbon and nitrogen dif- 
ferences as large as those in NGC 6752, the majority of 
giants in M53 have a smaller dispersion. 

4.2. Carbon abundances 

To calculate carbon abundances for the stars in our 
M53 sample, we compare the S{CH) values from the 
Kast spectra to those from the model spectra. Figure 7 
presents the data used for determining [C /Fe] , and shows 
S{CH) versus My for both data (circles) and models 
(solid lines). We interpolate linearly between the model 
isoabundance contours to calculate [C/Fe] for each star. 
The vertical error bars on the individual data points 
represent the observational uncertainty in S(CH). The 
mean value of that error is <7s(ch) = 0.007, which is 
equivalent to an error in [C/Fe] of roughly 0.05 dex. The 
derived values of [C/Fe] are accordingly rounded to the 
nearest multiple of 0.05 dex. Figure 11 shows [C/Fe] 
versus My for the CN-enhanced and CN-weak stars (as 
identified in Figure 3), and it shows a clearly declining 
trend, at least among the CN-weak giants, as is suggested 
by Figure 7. The carbon depletion rate in the full M53 
sample is 0.22 dex per V magnitude, which is similar 
to the rate measured in red giants in the halo field by 
Smith & Martell (2003). Studies of red giants in other 
globular clusters also show abundance trends that indi- 
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Fig. 10.— The 3883 A CN bandstrength index 5(3839) versus 
My for our sample of M53 red giants, with two sets of models over- 
plotted: the lower line is typical of a CN-weak giant with [C/Fe] set 
to —0.3 dex and [N/Fe] to 0.0 dex, and the upper line is typical of 
a CN-strong star, with [C/Fe]= -0.6 dex and [N/Fe]= +0.9 dex. 
The anticorrclated abundance difference between the two model 
tracks is taken from the Norris et al. (1981) and Da Costa & Cot- 
trell (1980) studies of NGC 6752, a typical CN-bimodal globular 
cluster. Some of the M53 stars fall on or near the model lines, but 
most lie between, indicating that they have intermediate carbon 
and nitrogen abundances. 



cate ongoing deep mixing (see, e.g., Carbon et al. 1982 
and Trefzger et al. 1983). The theoretical study of Denis- 
senkov & VandcnBcrg (2003), which models deep mixing 
as a diffusive process, finds that deep mixing rates among 
low-mass red giants ought to be fairly independent of 
mass and metallicity, which may explain why abundance 
trends are so similar among so many different groups of 
red giants. 

However, there is the possibility that the measured car- 
bon abundances are affected by a systematic error result- 
ing from initial model assumptions. Specifically, when 
the synthetic spectra were created, it was assumed that 
all stars had the same oxygen abundance regardless of 
luminosity. By not allowing for a decline in [O/Fe] with 
rising luminosity, our models pertain to a case in which 
deep mixing has not affected oxygen. 

Free oxygen in cool stellar atmospheres (T e g < 4500 
K) participates in molecular equilibrium by forming CO, 
which consumes carbon atoms that could otherwise form 
CH molecules. With a fixed carbon abundance, there- 
fore, an oxygen-enhanced star will have weaker CH bands 
than an oxygen-depleted star. If [O/Fe] starts at a value 
of +0.2 dex, but oxygen is depleted as stars evolve along 
the giant branch, due to some mixing process that our 
models do not take into account, then too large a carbon 
abundance will be inferred for the most evolved stars, 
since we will have assumed too large an oxygen abun- 
dance for them. This will act to make the decline in 
Figure 11 artificially shallow, and the true decrease in 
[C/Fe] with brightening My would be steeper than indi- 
cated. 

Briley et al. (1990) consider the effect of varying 
[O/Fe] from —0.6 to +0.4 dex on their [C/Fe] measure- 
ments for the low- metallicity cluster M55. Their Table 7 
shows that [C/Fe] determined from CH bandstrength can 
change by as much as 0.5 dex for a 1.0 dex change in as- 
sumed [O /Fe] , with the effect being largest in the coolest 
stars. Smith & Briley (2006) collected abundance data 
from the literature for the moderate-metallicity globular 
cluster M13. The data show a considerable spread in 
[O /Fe] at a given My , and a decline of roughly 1 dex in 
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Fig. 11. — Carbon abundance versus My for red giants in M53, 
with CN-enhanced stars (identified in Figure 3) plotted as filled 
circles. The decline in [C/Fe] with increasing luminosity expected 
from Figure 7 is seen here. As discussed in the text, it is possible 
that [C/Fe] is overestimated for the CN-strong stars brighter than 
My ~ —1.4 if overestimates of [O/Fe] are used in the models. 

the minimum [O/Fe] over the My range covered in our 
M53 sample. If oxygen behaves in an analogous way in 
M53, our use of constant- [O/Fe] models would cause sys- 
tematic errors in the [C/Fe] difference inferred between 
the faintest stars in our sample and those brighter gi- 
ants that have been diminished in oxygen. However, as 
noted in the previous paragraph, such an effect would not 
lead us to infer a spurious decline in carbon abundance 
along the red giant branch, but would act in the oppo- 
site sense to partially mask any carbon decrease that is 
present. Thus, the evidence of Figure 11 that evolution- 
ary changes do occur in the surface carbon abundances 
of the M53 red giants would remain valid, although the 
size of the effect may be underestimated in the present 
work. 

The hottest stars in our M53 sample should be the least 
affected by [O /Fe] variations because they are unlikely to 
have significant abundances of CO in their atmospheres. 
Suntzcff (1981) showed that relatively little CO forms 
above 4500 K in red giant photospheres at the metallic- 
ity of M13. According to the Victoria isochrone we have 
adopted for M53, a temperature of 4500 K corresponds 
to an absolute magnitude of My = —1.4. It is among 
stars fainter than this that the evidence for a CN-CH 
band anticorrelation in Figure 8 is greatest. If the at- 
mospheres of CN-enhanced stars have been subjected to 
greater CNO-processing than those of the CN-weak stars, 
they will have lower [O/Fe] and [C/Fe] abundances but 
higher [N/Fe] than CN-weak stars. This will produce 
a CN-CH band anticorrelation in stars hot enough to 
have little or no CO molecular formation in the atmo- 
sphere. However, in cooler CN-enhanced giants there 
would be less CO formation than in the CN-weak giants, 
which may partially counteract a CN-CH anticorrelation 
by leaving relatively more carbon available for CH for- 
mation. Figure 8, seen in this light, shows a CN-CH 
anticorrelation in the hotter stars that may be erased 
as CO formation sets in more strongly in the CN-weak 
stars. 

5. DISCUSSION 

The generalized histogram of 55(3839) shows that M53 
has a distribution of CN bandstrength whose width is 



greater than our 3-tr errors. We infer that M53 has a 
population of CN-enhanced giants, and not just CN-weak 
stars. However, there are two physical factors that erode 
the appearance of CN bimodality at low metallicity: first, 
in a cluster with a low [Fe/H], the low carbon and nitro- 
gen abundances can sustain only limited molecular CN 
formation. Second, deep mixing scenarios predict that 
mixing efficiency should be high at low metallicity (e.g., 
Sweigart & Mengel 1979; Charbonnel 1995), meaning 
that evolved stars in a low-metallicity cluster like M53 
will be depleted in carbon and enhanced in nitrogen from 
the dredge-up of C— >N and possibly O— >N-processed ma- 
terial. If the mixing continues to the point where carbon 
is less abundant than nitrogen, the abundance of CN 
molecules can decline with rising luminosity (see, e.g., 
the discussion of extremely C-depleted stars in Cohen et 
al. 2002 and Langer et al. 1985). 

The relatively small range in carbon and nitrogen 
abundance variation observed in M53 may also be an 
indication that the polluting material was not processed 
through the full CNO cycle. If this is the case, then 
M53 should not show an O-Na anticorrelation, because 
neither the O^N processing nor the NeNa-cycle process- 
ing that occurs at the same temperature will have taken 
place. The data collected for this study do not allow for 
an investigation into this question, but a followup study 
at higher resolution could be quite interesting. O-Na 
anticorrelations have been found in red giants in every 
globular cluster where they have been sought, meaning 
that their absence in M53 would be notable. 

In summary, we find a progressive decrease in [C/Fe] 
with brightening My among red giants in M53, indi- 
cating that deep mixing is occurring. The [C/Fe] and 
5(3839) measurements further reveal that some stars 
have carbon depletions and nitrogen enhancements (rel- 
ative to the cluster CN-weak stars) as large as those 
observed in the CN-strong giants of the classically CN- 
bimodal globular cluster NGC 6752 (although we do not 
see a strongly bimodal CN distribution in M53), indi- 
cating that whatever process is responsible for primor- 
dial abundance variations in NGC 6752 (as described in 
Caretta et al. 2005 and Suntzcff & Smith 1991) is also op- 
erating at the lower metallicity of M53 ([Fe/H]~ —1.8). 
We conclude that the rarity of CN-strong stars in other 
clusters with [Fe/H] ~ —1.8 or lower, such as M55, is 
more likely an effect of the overall low metallicity caus- 
ing diminished CN molecule formation and the reduced 
sensitivity of the index 5(3839) at low metallicity, and 
is not due to some unexplained lack of abundance inho- 
mogeneity 6 . A direct test of this hypothesis would be to 
obtain spectroscopy of the A3360 NH band for red giants 
in M53 and M55. In the case of M92, which is of even 
lower metallicity, it would seem that the formation of 
CN-strong giants is inhibited by the low overall metallic- 
ity of [Fe/H] 2.3 even though large [C/Fe] and [N/Fe] 

differences, as well as large [N/Fe] enhancements, were 
found among the red giants by Carbon et al. (1982). 

6 However, it is certainly the case that the fraction of CN-strong 
giants varies among globular clusters, as discussed in, e.g., Norris 
1987. The small number of M53 giants that lie near the upper 
isoabundance line in Figure 10 may indicate that M53 has a lower 
fraction of CN-enhanced stars than NGC 6752. 
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